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We review recent work on resonant Andreev tunneling through a strongly interacting quantum dot 
connected to a normal and to a superconducting lead. We derive a general expression for the current 
flowing in the structure and discuss the linear and non-linear transport in the nonperturbative 
regime. New effects associated to the Kondo resonance combined with the two-particle tunneling 
arise. The Kondo anomaly in the I — V characteristics depends on the relative size of the gap energy 
and the Kondo temperature. 



I. INTRODUCTION 

In recent years electrical transport through confined 
regions has seen impressive theoretical and experimental 
activity. Several well known phenomena have received 
renewed attention due to the present possibility of study- 
ing them in novel and more controllable situations. For 
instance, Andreev scattering (Andreev 1964), according 
to which a particle-like excitation impinging on a nor- 
mal metal-superconductor interface is reflected back as a 
hole-like excitation, has been shown to be the key mech- 
anism controlling transport in a variety of hybrid meso- 
scopic superconducting devices. (For a review see, for 
example, Hekking, Schon, and Averin 1994, Beenakker 
1995, Lambert and Raimondi 1998). At the same time, 
electron-electron interaction in small confined regions, or 
quantum dots (QDs), has been shown to lead to the so- 
called Coulomb blockade of electrical transport. This 
occurs at low temperatures, when the Fermi energy of 
the contacts falls in the gap between the ground state 
energies of the dot with N and N + 1 electrons (Kastner 
1992). However, a QD attached to metallic leads resem- 
bles an impurity level in a metal. As a consequence, even 
in the Coulomb blockade regime, transport will occur due 
to the Kondo effect (Glazman and Raikh 1988, Ng 1988). 
This is due to the formation of a spin singlet between the 
impurity level and the conduction electrons, which gives 
rise to a quasiparticle peak at the Fermi energy in the 
dot spectral function. This suggestion has been explored 
theoretically by several groups (Meir, Wingreen, and 
Lee 1991,1993, Ng 1993, Levy Yeyati, Martin-Rodero, 
and Flores 1993, Schoeller and Schon 1994, Hettler and 
Schoeller 1995). This has lead to the prediction of a 
zero-bias anomaly in the current voltage characteristics 
and an increase of the linear conductance in the Coulomb 
blockade regime for decreasing temperature. Such phe- 
nomena have indeed been observed in different QD sys- 
tems (Ralph and Buhrman 1994, Goldhaber-Gordon et 
al. 1998, Cronenwett et al. 1998). 

What will happen if the QD is coupled to a normal and 
a superconducting lead as shown schematically in Fig^? 
Does the zero-bias anomaly observed in the N-QD-N case 



survive in the N-QD-S case? Such a problem has been 
investigated recently by various groups (Fazio and Rai- 
mondi 1998, K. Kang 1999, Schwab and Raimondi 1999, 
Sun, Wang, and Lin 1999, Clerk, Ambegaokar, Hersh- 
field 1999). In this paper we review mainly the work 
done in Refs. (Fazio and Raimondi 1998) and (Schwab 
and Raimondi 1999). 

The problem one is faced with has all the difficulty of 
the original Kondo problem. It is well known that in 
the limits of high and low temperatures (here high and 
low are with respect to the Kondo temperature), qualita- 
tively correct results are obtained by means of different 
techniques. Such an attitude we adopt here, having in 
mind to analyze the interplay between Andreev scatter- 
ing and Kondo physics. We use a simple equation-of- 
motion approach in the perturbative regime above the 
Kondo temperature. This corresponds to taking into 
account only the leading logarithmic corrections in the 
renormalization group sense. The cross-over behaviour 
to low temperature may be rather complicated, but as 
happens in the standard Kondo problem, we expect that 
the extreme low-temperature phase can be described in 
simple terms. Such a description is obtained by means of 
the slave-boson technique in a mean-field approximation. 
In the high temperature regime, we find a suppressed An- 
dreev current at low bias voltage due to the competition 
between the Coulomb energy and the superconducting 
proximity effect in the QD. In the low-temperature case, 
in contrast, our analysis predicts that the linear conduc- 
tance of the N-QD-S system is enhanced as compared to 
the normal case and may reach the maximum universal 
value of Gns = 4e 2 //i which is twice the maximum for 
the N-QD-N system. 

Before entering the technical details, it is also worth- 
while discussing the relevant energy scales. In investigat- 
ing Andreev scattering at a normal metal-superconductor 
interface one is often interested in voltages and temper- 
atures well below the energy scale set by the supercon- 
ducting gap A and one could be tempted to take A — > oo 
from the outset. However, the dot charging energy U (see 
below) introduces another energy scale into the problem, 
and the physics is different in the two limits U A 
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and A 3> U. In order to appreciate this point, let us 
consider the Hamiltonian of an isolated QD, plus a term 
which models Andreev scattering in the limit A — > oo 
by Ha — t A d^d^ + c.c. For this Hamiltonian the off- 
diagonal element of the inverse dot Green's function reads 
-t A {l + U(e + e d + U)/[e 2 -(e d + U) 2 ~t A }} and vanishes 
in the U — > oo limit, i.e., in the limiting sequence A — » oo 
first, U — > oo after, the induced superconductivity and 
hence transport in the dot is completely suppressed. In 
order to have Andreev scattering in the large U limit two 
electrons have to enter the superconductor without dou- 
bly occupying the QD. This can only happen on a time 
scale of the order 1/A. We therefore concentrate on the 
limit U ^> A in the following analysis. 

Our paper is organized as follows. In Section II we 
derive a formula for the current through the dot in the 
presence of Andreev scattering. Calculating the current 
explicitly requires the dot Green's function. In Sections 
III and IV we present the equation-of-motion and slave- 
boson approaches, respectively. A few conclusions are 
drawn in section V. 

II. THE ANDREEV CURRENT FORMULA 

In this section we define the system under consider- 
ation and derive an expression for the Andreev current 
which holds in the presence of electron-electron interac- 
tion. 

The model Hamiltonian for the N-QD-S system is 

H = H N + H g + H D + i?x,N + ^t,s (1) 

where 7?n, Hs, are the Hamiltonians of the normal and 
the superconducting leads (A being the superconducting 
gap), 

k,cr 

H S =Y1 e k4 k(T c S!k(T + X!( A 4kT c s,-ki + c - c -) • ( 3 ) 

k,CT k 

If we restrict ourselves to temperatures and bias voltages 
much smaller than the average level spacing in the QD, 
then transport occurs through a single level. In this case 
the Hamiltonian of the quantum dot Hd reads 

H D = e d d\d a + Un d ^n dl . (4) 

The level e d is assumed to be spin degenerate and the 
electron-electron interaction is included through the on- 
site repulsion U (~ 1 — 5K for currently fabricated QDs). 
Experimentally the position of the dot level can be mod- 
ulated by an external gate voltage. Tunneling between 
the leads and the dot is described by i?x,N an d -Ht,s 



kcr 

where r\ = N, S and is the tunneling amplitude. 

The starting point for deriving the current formula is 
I = — e(d/dt)(N^) = ie [JVn, H] where is the electron 
number operator in the normal lead. In the case of a 
hybrid structure, like the one we are considering now, it 
turns out to be more convenient to evaluate the current 
in the normal lead. The average current can be rewritten 
as follows 

7 = 2eIm^V N (4 kCT rf ff ) . (6) 

Since we deal with a nonequilibrium situation we work in 
the framework of the Keldysh technique, as employed in 
the literature (Meir and Wingreen 1992). By introducing 
the Nambu notation 

the average current in eq.(^) requires the evalua- 
tion of the lesser Green's function G^g k (t,t') — 

'(^pk^l^W) which, by means of the Dyson equa- 
tion, can be expressed in terms of the exact Green's func- 
tion G of the QD and the free Green's function of the 
normal lead (?N.k- The lesser component of the Dyson 
equation reads 

+G<(e)^ N ^ W k (e) . (8) 

where G R ( A \ G< are the retarded (advanced) and 
the lesser Green's functions of the dot (for example 
G R {t) = — i6»(t)<{^(t), 0+(O)}>). Using the relation G< = 
G R T, < G A and the expression for the Green's functions 
in the normal lead (diagonal in Nambu space) 5n k ( e ) = 
l/(c -e k + i V ), <?< k (e) = 27ri/(e)5(e - e k ) (f(x) is the 
Fermi function), it is possible to rewrite the current in 
the following form 

1 = ie f-oo #r N Tr{f z G*(e)[X:«(e)/ N (e) - 

h(e)E A (e) + t<(e)]G A (e)} (9) 

In the formula above, we introduced the elastic rate 
Tn(c) = 27ry^ k |Vn| 2 £(£ — e k ). The diagonal matrix /n 
has elements / N ,n = f(e - /x N ) and / N ^2 = 1 - f(-e - 
/Un)- The normal electrode is kept at a chemical poten- 
tial /iN = —eV while that of the superconductor is fixed 
to zero. 

Up to now no approximations were involved. In or- 
der to determine the current an expression for the self- 
energy of the QD should be found. A determination of 
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I requires both the lesser and retarded parts of the self- 
energy. We formulate an ansatz for the lesser Green's 
function which is expressed solely in terms of the retarded 
one. This ansatz automatically guarantees current con- 
servation. This generalizes an ansatz put forward by (Ng 
1996) for two normal leads to the calculation of S < (e) in 
the presence of a superconducting lead. 

In order to understand how the ansatz is constructed 
it is useful to start from the noninteracting dot. In this 
case, £ <:> (e) can be computed exactly and it is expressed 
in terms of the retarded and advanced self-energies as 



5£(e) 



E 

N.S 



E£„(e)/„(e) 



/,(e)S^(e) (10) 



y> 



(e) = -E J7 =N,s{^(e)[i-/,(e)] 
[i-/ 7) (e)]E^(e)}. 



(11) 



In this case (with no interaction) the nonequilibrium self- 
energy has the same form as in equilibrium but with the 
Fermi functions of the two leads kept at different chem- 
ical potentials. The idea is to assume that even in the 
presence of interaction the dependence on the Fermi dis- 
tribution is the same and that both the lesser and the 
greater functions depend on a single function A such that 



1995) for the non-interacting case to the case of a strongly 
interacting dot . 

In the following sections we apply this formula to the 
study of the I — V curves of a N-QD-S system. 



III. EQUATION OF MOTION 

Our aim is to find an equation for the dot's Green's 
function G{t) = -i(T [</>{t) 0(0)]). The retarded and ad- 
vanced Green's functions necessary for the evaluation of 
the current can easily be obtained from the time-ordered 
one. In this section we will describe the equation-of- 
motion approach. It is useful to introduce quasiparticle 
operators by performing a Bogolubov transformation 



7 V M = "^kc^ko- + sign(cr)e lip v v ^c f 



T}, — k(T 



7^,-kF = -sign(cr)e I< <X k c, ):k(T + u vM cl h _ k - 



(16) 
(17) 



w here u^ k (v ^ k ) = (1 ± e k/ %, k ) /2 , E nM = 
\/e k + I \ | 2 , and A r) = e i(p \ A v \. In the normal lead 
we have An = and the Bogolubov transformation re- 
duces to the identity transformation. The hopping term 
becomes in the quasiparticle basis 



E<A 



£>A 



The matrix A is determined by the condition 



£ < — £ > = £ 



R 



(12) 



(13) 



which is a general property of the Keldysh Green's func- 
tions. As already mentioned, this ansatz leads to cur- 
rent conservation. Moreover it is exact both in the non- 
interacting limit, U = 0, and in absence of superconduc- 
tivity, A = 0. As a result 



A\-l, 



(14) 



Eq. ( |l4|) allows us to evaluate the expression of the current 
(^) , once we know the retarded Green's function of the 
dot. The expression for the current can be greatly sim- 
plified in the relevant limit U, A ^> fcgT, eV. In this case 
the non-interacting self-energy due to the superconduct- 
ing lead E^g" 4 ' is real and purely off-diagonal whereas 

that due to the normal lead, E^n, is diagonal. Us- 
ing these properties of the self-energy and substituting 
eq.(|lj) in eq.(||) we obtain the following form for the 
Andrcev current through a QD. 



I = ie 



£^r N Tr (r,G«(e) [s«(e), / N (e)] G A (e)} 



(15) 



This equation generalizes the formula obtained by 
(Beenakker 1992, Claughton, Leadbeater, and Lambert 



T r)i k — 



?7,k 



(18) 



The Nambu formalism is introduced for the 7 operator 
via ^.k = (7»j,kTi7^ - k j.)- We start by writing down the 
equation of motion for the operators <f> and ^: 



r;,k 



d ~ - 



(19) 
(20) 



In the equation for (f>, the new operator $ = 
(dld^d^, djdjdf), appears because of the interaction. It 
is straightforward to iterate once the equation of motion 
and get an equation for $ 



(21) 



where 



N = 



d\d], did 



d\d\, 



dtdf 



Here the superscript t means the transpose. In eq.(|2l|) 
again several new operators have appeared. In order to 
get a closed system of equations one has to truncate the 
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hierarchy at some point. Following the decoupling proce- 
dure used in the absence of superconducting leads (Meir, 
Wingreen and Lee 1991), we neglect correlations in the 
leads. To see how this is achieved, consider the general 
structure of all the new operators entering eq.(pl|) for $. 
They are products of two d and one 7 operator. By iter- 
ating once more the equation of motion one obtains new 
terms which contain two 7 and one d operator. The ap- 
proximation consists then in replacing a pair of two such 
7 operators with their statistical average. In such a way 
all the new operators entering the eq.(|2l]) for $ can be 
expressed in terms of $ itself and <f>. A further approxi- 
mation is done by considering the limit of large Coulomb 
interaction, i.e., U — > 00. In this case one can safely 
neglect, in the third order of iteration of the equation of 
motion, all operators having two creation or two annihi- 
lation d operators, because they will give rise to terms of 
the order 1/U. We write here the resulting equations of 
motion for G and G n = -i (T [$(t)0(O)]). They are 

e d 6- 2 ) G(uj) = a + t a {Lu)G{Lu) + Ua z G H (cu) (22) 



(ujfr 



Ufr z G u (uj) = -(N) + E/(w)G(u>) 



(23) 



and the equation for the single particle Green's function 
becomes 



(cjo-q - e d fr z - fr(u>)) G(uj) =fr - (N), 



(24) 



where a = £0 + The noninteracting self-energy £0 is 
obtained as 



A 2 

v 



lu A n 
-A v -oj 



(25) 



/ w 

with r j; = 2tt J^k I 1 2 <5(e — ek). For the calculation 
of the interacting self-energy £/ (uj) we assume that in 



the superconducting reservoir there are no quasiparticles 
present (we consider energies much smaller than the gap 
A) (7s ko-Ts.ktr) = 0, and that in the normal reservoir 
the quasiparticle population follows the Fermi distribu- 
tion (7J, ko .7N,k ( r} = / (ek - Mn)- We then have 



"+-Ek.s 



/(<* 



2Mk,S"k,S 

"--Ek,s 
"k.sl 2 







(26) 



(27) 



In what follows, we consider W ^> A ^ T, u> with W 
the bandwidth. The self-energy Si,s,ii is weakly en- 
ergy dependent since u> A < £"k.s- The imaginary 
part of the diagonal elements of the self-energy Ss van- 
ishes, so quasiparticles present in the dot cannot decay 
by tunneling into the superconductor. As a result the 
contribution to the self-energy due to the superconduct- 
ing lead simply shifts the dot level to the new value 
(d ~ (d + (rs/27r) In W/A. The divergence of the level 
energy renormalization with A reveals that the process 
occurs via a virtual state in which a quasiparticle is cre- 
ated in the superconductor. Substituting the expression 
of the QD's Green's function in eq.(|l5|) we get the result 



I(V) = 



f 

J — ( 



de 



f(e-eV)-f(e + eV) 
2e 



G NS (e) (28) 



with 



G NS (e) 



4e 2 



2(r N r s ( e )) 2 [(r liN + r 2 , N )/2r N )] 



h [4(e-ei)(e + e 2 ) -ri, N r 2 , N 



r 2 (e)] 2 + 4 [r 1;N (e + e 2 ) + r 2 , N (e - £i)] 2 



(29) 



where e 1(2) = e d ± Re<7ii( 22) , r 1(2) , N = -2Im<7 U ( 2 2), 
and T s (e) = 2Re<7i 2 = r 5 2e/(7rA). Notice that the 
anomalous non- interacting £0,12 an( i interacting £1,12 
self-energies exactly cancel in the zero energy limit, re- 
sulting in a linear behaviour for energies smaller than the 
gap Q. The spectral function Gns(e) associated with 
the resonant Andreev tunneling is plotted in Fig. |^ for 
various bias voltages. Several features are worth noticing. 
First, two peaks at ±e are due to particle and hole bare 
levels. Note that in the interacting case the bare level 
energy includes the renormalization due to the supercon- 
ducting electrode self-energy as discussed above. Second, 
at low temperatures the spectral function at the Fermi 
energy is completely suppressed, in contrast to what hap- 
pens in N-QD-N case. Quite remarkably, at finite pos- 
itive (negative) voltages a Kondo peak develops pinned 
to the Fermi level of the normal metal while a small kink 



develops at negative (positive) voltages. At finite volt- 
ages hole and particle energies differ by 2eV, and while 
the electron (hole) is on resonance for positive (negative) 
voltage, the Andreev reflected hole (electron) is off reso- 
nance with respect to the shifted Fermi level. 

The differential conductance, for various temperatures, 
is shown in Fig. [|. Lowering the temperature a zero- 
bias anomaly starts to develop where the conductance 
is strongly suppressed at low voltages. From cqs.(p8|- 
|29| ) we conclude that the linear conductance is roughly 
proportional to T 2 /A 2 , so that it seems to be com- 
pletly suppressed in the zero temperature limit. However, 
the equation-of-motion approach is quantitatively reli- 
able only above the Kondo temperature. Hence the anal- 
ysis carried out here applies to the regime A > T > Tk- 
An interesting question to ask is what happens in the op- 
posite regime when the Kondo temperature dominates, 
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i.e., when Tk > A > T. This regime cannot be explored 
by the equation-of-motion approach. For this reason in 
the next section we will investigate the low temperature 
regime by means of the slave-boson technique. 



b 2 + Tj2^[0(^ n )a z 
b X + bpTy^Tr ^(ie„)f (ie„) 



= 0, 
= 0, 



(32) 
(33) 



IV. SLAVE BOSON MEAN FIELD 
APPROXIMATION 

In this section we extend the analysis of the previ- 
ous section to the extreme low temperature regime. To 
this end we use slave-boson mean field theory. This ap- 
proach has been successfully applied to the low temper- 
ature properties of a Kondo impurity in the presence of 
normal (Barnes 1976, Coleman 1984, Read and Newns 
1983, Read 1985) as well as for superconducting conduc- 
tion electrons (Borkowski and Hirschfeld 1994). Despite 
its simplicity, this method captures the main physical as- 
pects of the Fermi liquid regime at low temperatures, i.e., 
the formation of a many-body resonance at the Fermi en- 
ergy. For this reason it presents a convenient framework 
in which to study the interplay between Andreev scatter- 
ing and Coulomb interactions. 

Again we consider an infinite on-site repulsion U, so 
processes where the dot level is doubly occupied are ex- 
cluded. The dot level is represented as d\ — where 
the fermion f a and the boson b describe the singly occu- 
pied and empty dot states. Since the dot is either empty 
or singly occupied, the constraint + flf a = 1 has 
to be fulfilled. 

In mean field approximation, the operator b is replaced 
by a c- number bo, and the constraint is fulfilled only on 
average. This is achieved by introducing a chemical po- 
tential Ao for the pseudo particles. Notice that one ends 
up with a non-interacting-like problem with renormal- 
ized parameters, i.e., an energy shift for the dot level 
(d —>■ £d + Ao = e<j and a multiplicatively renormalized 
tunneling amplitude Vq — > boV n . 

We discuss the mean field equations and their so- 
lution first in equilibrium and then generalize to non- 
equilibrium. We start from the impurity part of the free 
energy, which in the presence of both normal and super- 
conducting leads is given by 

F = -T^Trm[ie„<7o - e d a z - 6^r(ie„)] + \ bl + e d - M, 

(30) 

where e„ is a fermionic Matsubara frequency, <j l are the 
Pauli matrices, and 



r(ie n ) = ^ \y^o z g ri .yX\(. n )a z 



(31) 



with g^.k being the Green's function of the lead r\. 

By minimizing the free energy with respect to Ao and 
bo we find the equations 



which have to be solved self-consistently. G(ie n ) is 
the pseudo fermion Green's function given by G(ie n ) — 
[ie n ao — £d<J z — b^T (ie n )] . Both in the limit of small 
and large superconducting gap, we are able to solve the 
mean field equations analytically as demonstrated here 
below. The first equation, eq. (|32]) , is the constraint. 
Since the pseudo fermion level is at maximum singly oc- 
cupied, the renormalized level is above the Fermi energy. 
In the Kondo limit, where the occupancy is nearly one, 
we find that < < b^T^ + Ts), i.e. A « |e<j| and 
£d ~ 0. The renormalization of the tunneling amplitude 
is determined from eq. ([33]) . A trivial solution bo = 
always exists. The solutions which minimize the free en- 
ergy, however, are those with bo ^ 0. By introducing a 
flat density of states in the leads and the tunneling rates 
r j( = 2irNo\V v \ 2 , the elements of the matrix r(ie„) are 
Tn = T 2 2 = -171 and r 12 = = 72, where 



71 = signe- 



rs 



A 



2 y^T^p 



72 



2 y/e* + |AP 
(34) 



Restricting ourselves to zero temperature, we replace the 
Matsubara sum in eq. ( [33] ) by an integral and obtain 



kdl = 4 



w de 7l ( e + hg 7l ) + frg| 72 F 

'0 27T ( e + &2 7l) 2 +6 4| 72 | 2 



(35) 



where IF is a cut-off of order the band-width. We sim- 
plify the integral by approximating 71 and 72 as 



r N /2 for e < A 

~ ' (r N + r s )/2 for e > A 



72 



r s /2 for e < A 
for e > A 

(36) 



The result is 



£n (2A 
2n ln 



b A Y 2 
u a L s 



r N + r £ 

2vr 



■ In ■ 



4IF 2 



(2A + b 2 r s + b 2 r N y 



(37) 



where we neglect a term proportional to Ts , but without 
any logarithmic factor. If A is much smaller than the 
Kondo temperature which is given by Tk = ^o^n + ^o^s > 
A is negligible. One can then easily solve eq.(^) for 
b^ and obtain the result for two normal leads with total 
tunneling rate Tn + Ts : 



foo(r N + r s ) = 2FFexp 



(38) 
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In the opposite limit, where A is much larger than Tk, 
we find 

(39) 



The results agree qualitatively with what we expect from 
scaling arguments for the Anderson model. In the per- 
turbative regime, a logarithmic correction to has been 
found (Haldane 1978). This applies to the case of a large 
gap, since scaling due to the superconducting electrons 
stops at energies of the order A, giving rise to a log- 
arithmic renormalization of ed, as seen in eq. (|39|) . In 
the case of a small gap, the superconducting lead con- 
tributes to scaling down to low energies, where one en- 
ters the strong coupling regime. Presumably, the fixed 
point is still reached for energies of the order of Tk, much 
greater than A, so that the Kondo temperature does not 
depend on A, as indeed found in cq.(|38|). Notice that 
in the presence of normal electrons, we always find a 
non-trivial solution of the mean field equations. This is 
to be contrasted with the case of superconducting elec- 
trons only, Tn = 0, where for large gap only the solution 
bo = exists, and there is no Kondo effect (Borkowski 
and Hirschfeld 1994). 

In a non-equilibrium situation, when a voltage is ap- 
plied between the two leads, the mean field parameters 
cannot be obtained by minimizing the free energy. How- 
ever the mean field equations (^J^t ) can also be derived 
using a self-consistent diagrammatic method (Millis and 
Lee 1987). Then it is straightforward to generalize to 
non-equilibrium. The equations read 



, de 

11 



Q < (e)a z =0 (40) 



A &o - ifro 



de 

2n~ 



Tr 



g R (e)f<(e) + g<(e)f A (e) =0, (41) 



where the lesser Green's function Q < (t,t') = {t')(f)(t)) 
has been introduced, with cf> = (/f,/f). The lesser and 

advanced matrix T is defined in analogy to its equilib- 
rium version in eq.(|3l|). To obtain Q < , we use the gen- 
eral relation Q K — Q R T, < Q A , where at mean- field level 



b 2 n T< and 



?7,k 



V r . 



(42) 



Note that the superconducting lead does not contribute 
to £<(e) for |e| < A. 

We have solved the mean-field equations in the pres- 
ence of an external voltage numerically. As long as 
\eV\ <C Tk the solution is almost independent of the 
voltage. For large voltage, \eV\ 3> Ik, we have found 



that the Kondo peak is pinned to the chemical potential 
in the normal lead, i.e. ed —> e^ — eV, and the peak width 
is decreased. 

The Andreev current can now be determined using the 
current formula eq. (|l5|) , with the dot Green's function 
G - 

eq.(|2S|), is determined as 



b^Q. The spectral function, which is defined as in 



G NS (e) = 



4e 2 



4(r N r s ) 2 



h (4e 2 - 4e 2 "In - TlY + 161* e? 



(43) 



Here the renormalized tunneling rates Ts,n = ^o^n, 
and e = e(l + 6qTs/2A) are introduced. One recovers 
the current formula for a non-interacting quantum dot 
(Beenakker 1992), with renormalized parameters which 
are voltage dependent. On resonance, when id ~ and 
e = 0, the small renormalization factor bo drops out. The 
differential conductance becomes maximal when Tn = Tg 
with Gns, max = Ae 2 /h, twice the maximum for a N-QD- 
N system. For large voltage Gns drops quickly, since 
the resonance moves away from zero energy, |ed| ~ |eV|. 
Gns as a function of energy is proportional to (eV)~ 2 
near e « ±ey and proportional to (eV)~ A near e = 0. As 
a consequence the current decreases with increasing volt- 
age, leading to a negative differential conductance. This 
is also demonstrated in Fig. ^, where the current as a 
function of voltage obtained using Gns of eq.([43"|) with the 
numerically determined, voltage dependent mean field 
parameters. The results were obtained with e d = W/3, 
T N = T s = 0.14W, and A = 0.2W. For low voltages the 
current is, to good approximation, given by I — Ae 2 V/h, 
whereas the current drops when the voltage exceeds the 
Kondo temperature. 

Finally, we want to comment on the reliability of our 
results. The success of slave-boson mean field theory 
stems from the fact that it captures the Fermi-liquid 
regime at low temperature. If the N-QD-S system scales 
to a Fermi liquid at low temperature, Gns as given in 
eq.(^) is exact in the low temperature, low voltage limit. 
Since it is known that slave boson mean field theory has 
problems in describing dynamical properties, the results 
far away from equilibrium need to be treated with cau- 
tion. 

Within the Fermi-liquid point of view, the present 
mean field approach allows us to estimate the parameters 
entering eq.(^3|). In particular, we found that Tn and Tg 
renormalize equally, but this may no longer be the case 
when considering higher order corrections. For illustra- 
tion, we estimate the effect of residual quasiparticlc in- 
teraction in the limit A <C Tk- By assuming an effective 
quasiparticle interaction of the form H- lrA — Un^ni, we 
find to first order in U no corrections to Tn, while, as 
one could have expected, repulsive quasiparticle interac- 
tion suppresses the renormalized coupling to the super- 
conductor, fg = 6^r s [l - 4(C//7rr K )(A/TK)lnT K /2A]. 



G 



V. CONCLUSIONS 

We have studied Andreev tunneling in a normal metal- 
quantum dot-superconductor device, and obtained a gen- 
eral formula for the current through the device. We 
calculated the current voltage characteristics within the 
equation-of-motion approach, where we found a suppres- 
sion of the conductance at low temperature. However it 
is known that this approach does not provide quantita- 
tive results in the Kondo regime. Therefore we extended 
the analysis to the extreme low temperature regime using 
slave boson mean field theory. In the regime when the 
superconducting gap is smaller than the Kondo tempera- 
ture, we found an enhanced Andreev current at low bias 
voltage due to the Kondo effect. The zero bias conduc- 
tance is maximum with the universal value Gns = 4e 2 //i 
when the renormalized tunneling rates Tn and Ts are 
equal. We identified the ratio A/Tk as an important 
parameter. In the case A <C 7k, the Kondo resonance 
forms as for two normal leads. The condition Tn = Tg 
coincides with equal bare tunneling rates Tn = Tg . In the 
case of large gap, quasiparticle interaction suppresses Fs, 
nevertheless the conductance maximum condition may be 
achieved in an asymmetric QD where Ts > Tn- 
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FIG. 1. The system under consideration. A quantum dot 
coupled by tunnel barriers to a normal and to a supercon- 
ducting electrode. The position of the level in the dots can 
be tuned by means of the gate voltage V g 
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FIG. 3. The differential conductance of the N-QD-S de- 
vice, in units of 4e 2 /h, is plotted for different temperatures 
(T = 0.0001 solid line, T = 0.001 dotted line, T = 0.01 
dot-dashed line, e~ d = -0.04, T s = T N = 0.02, A = 0.1, in 
units of the bandwidth W). 
0.008 
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FIG. 2. The spectral density for the two particle tunneling 
Gns(e) is plotted for various bias voltages ( V = solid line, 
V = -0.015 dotted line, V = 0.015 dashed line, e d = -0.07, 
r s = T N = 0.02 and T = 0.0001, A = 0.1, in units of the 
bandwidth W). In the inset the same curves are shown in an 
extended scale. 
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FIG. 4. Current voltage characteristics of a quantum dot 
at zero temperature as obtained within slave boson mean field 
theory. e d = 1/3, r N = Ts = 0.14, and A = 0.2 in units of 
the bandwidth. The Kondo temperature in equilibrium is 
T K ~ &or N = 0.014. 
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